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Abstract. A quantitative analysis of time series of ultraviolet 
spectra from a sample of 10 bright O-type stars (cf. Kaper et al. 



1996, Paper I) is presented. Migrating discrete absorption com- 



ponents (DACs), responsible for the observed variability in the 
UV resonance doublets, are modeled. To isolate the DACs from 
the underlying P Cygni lines, a method is developed to con- 
struct a template ("least-absorption") spectrum for each star. 
The central velocity, central optical depth, width, and column 
density of each pair of DACs is measured and studied as a func- 
tion of time. 

It turns out that the column density of a DAC first increases 
and subsequently decreases with time when the component is 
approaching its asymptotic velocity. Sometimes a DAC van- 
ishes before this velocity is reached. In some cases the asymp- 
totic DAC velocity systematically differs from event to event. 

In order to determine the characteristic timescale(s) of DAC 
variability, Fourier (CLEAN) analyses have been performed on 
the time series. The recurrence timescale of DACs is derived 
for most targets, and consistent results are obtained for differ- 
ent spectral lines. The DAC recurrence timescale is interpreted 
as an integer fraction of the stellar rotation period. In some 
datasets the variability in the blue edge of the P Cygni lines 
exhibits a longer period than the DAC variability. This might 
be related to the systematic difference in asymptotic velocity 
of successive DACs. 

The phase information provided by the Fourier analysis 
confirms the expected change in phase with increasing velocity. 
This supports the interpretation that the DACs are responsible 
for the detected periodicity. The phase diagram for the O gi- 
ant £ Per shows clear evidence for so-called "phase bowing", 
which is an observational indication for the presence of curved 
wind structures like corotating interaction regions in the stellar 
wind. An important difference with the results obtained for the 
B supergiant HD 64760 (Fullerton et al. |1997| ) is that in this 
O star the phase bowing can be associated with the DACs. No 



other O stars in our sample convincingly show phase bowing, 
but this could be simply due to the absence of periodic signal 
and hence coherent phase behaviour at low wind velocities. 

Key words: Stars: early type - Stars: magnetic fields - Stars: 
mass loss - Stars: oscillations - Ultraviolet: stars 



1. Introduction 

Variability is a fundamental property of the radiation-driven 
winds of early-type stars. Discrete absorption components 
(DACs) are the most prominent features of wind variability. 
They migrate from red to blue in UV P Cygni lines, narrowing 
in width when they approach the terminal velocity of the 
stellar wind. Obviously, DACs cannot be observed in saturated 
P Cygni profiles; however, the steep blue edges of these pro- 
files often show regular shifts of up to 10% in velocity. Edge 
variability is probably related to the DAC behaviour, but the 
precise phase relation has not been unraveled yet. For a more 
extensive introduction and various examples of the DAC phe- 
nomenon we refer to the first paper in this series (Kaper et al. 



1996, Paper I). Reviews on this subject were presented by Hen- 
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richs ( |1984|1988| ), HowarthJ199J), Kaper & Henrichs (|1994j), 
Prinja ( |l998| )Tand Kaper ( |1998| ). 

One of the key problems is to understand why DACs start 
to develop and how they evolve. Since DAC behaviour is dif- 
ferent for different stars, and even shows detailed changes from 
year to year for a given star, a detailed quantitative description 
of the observed variability over a significant amount of time 
(years) is essential. It should be emphasized that in the plan- 
ning phase of this project the sampling times of the targets have 
been carefully tuned to each individual star, based on sample 
spectra collected in earlier pilot studies. As a result, in Paper I, 
time series of more than 600 high-resolution ultraviolet spectra 
of 10 O-type stars were presented. In this follow-up paper we 
analyse the series of spectra obtained with the International Ul- 
traviolet Explorer (IUE) in a quantitative manner and evaluate 
the behaviour of DACs and edge variability in the UV P Cygni 
lines. In the following section we describe the data analysis and 
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Table 1. Main observational properties of the program stars. Notes: (a) Walborn 1972, except HD 210839 Walborn 1973; (b) 
Penny 1996; (c) Prinja et al. 1990 A typical uncertainty in is 50 km s -1 . For comparison, also the (maximum) asymptotic 



velocity of the DACs, w asym p is given based on our datasets. In the last column we list the number of highest flux points selected 
from the total set of spectra to construct the template spectrum. 
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present our methods to isolate and model the DACs in the ultra- 
violet spectra. Furthermore, a period-search analysis was per- 
formed to derive the characteristic timescales of both DAC and 
edge variability. The results of the modeling and the time se- 
ries analyses are presented in section 3. We compare the results 
obtained for different stars in section 4 and draw general con- 
clusions. In the last section we discuss the observed variability 
and its regular behaviour in terms of the corotating interaction 
regions model of Cranmer & Owocki (1996). 



2. Data analysis 

High-resolution (R =10,000) ultraviolet spectra of the 10 stud- 
ied O-type stars (see Table 1 for our sample) were obtained 
with the Short Wavelength Prime (SWP) camera on board the 
IUE satellite. Time series of the relevant UV P Cygni lines are 
presented in Paper I, together with a detailed description of 
the targets and their observational history. For each of the 10 
O stars, the high-resolution SWP spectra form a homogeneous 
dataset resulting from similar observational constraints and a 
uniform reduction procedure (cf. Paper I). 

2.1. Construction of template spectrum 

To isolate the migrating DACs from the underlying P Cygni 
profiles, a template (in this case a "least-absorption") spectrum 
is needed that represents the undisturbed-wind profile to which 
the P Cygni profile apparently returns after the passage of a 
DAC (e.g. Prinja et al. |1987 , Paper I). Subsequent division of 
the individual spectra by this template then results in quotient 
spectra, which are used to model the DACs (see next subsec- 
tion). 

Close inspection of the intrinsic variability in the P Cygni 
lines (cf. Paper I) suggests that the observed changes are mainly 
due to variations in the wind column in front of the stellar disk, 
since the wind variability is only found in the blue-shifted ab- 



sorption troughs, and not in the P Cygni emission peaks. On 
the one hand, this is what one would expect in the case of rela- 
tively modest variations in the density and/or velocity structure 
of the stellar-wind. To illustrate this, one could imagine a small 
blob of gas in the stellar wind absorbing photons at a particular 
velocity and reemitting them isotropically. When this blob is in 
the line of sight, it leaves a blue-shifted absorption feature, so 
that we can "count" all the photons that were scattered out of 
the absorbing column. For a blob outside the line of sight, only 
the few photons that are reemitted in the direction of the ob- 
server can be detected. On the other hand, the absorbing blob 
has to cover a significant fraction of the stellar disk in order to 
give rise to an observable change in the P Cygni absorption. 

Our basic assumption is that the changes in the P Cygni 
profiles are caused by variable amounts of additional absorp- 
tion caused by material in the line of sight. In terms of Sobolev 
optical depth, a change in P Cygni absorption could be due to a 
(local) change in the wind density or the velocity gradient. The 
simplest construction of a reference template spectrum is to se- 
lect the highest flux point (per wavelength bin) from the avail- 
able spectra. Such a procedure applied to non-variable parts 
of the spectrum would, however, yield a template that is sys- 
tematically too high because of instrumental and photon noise. 
Here we describe a method which automatically corrects for 
this overestimation. 

Our method is based on a paper by Peat & Pemberton 



( 197C ) which documents computer programs for the automatic 



reduction of stellar spectrograms, in particular to locate the 
continuum level in a spectrogram. The large number of avail- 
able spectra (N, listed in Table 1) enables the construction of a 
"noise-corrected" template spectrum. We assume that the flux 
points not disturbed by DACs are normally distributed with 
mean value /iq(X) and variance (J^(X) within each wavelength 
bin (taken 0. 1 A) at wavelength A. We choose a suitable num- 
ber n (last column in Table 1), not too small, but small enough 
to have the n highest flux values in every bin sufficiently free 
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Fig. 1. To demonstrate our method for the construction of a 
"least-absorption" template, we show the Si IV resonance dou- 
blet of the 07.5 giant 68 Cyg. The dotted line indicates the av- 
erage /i (upper panel) and standard deviation a (lower panel) 
of the 5 highest flux points (out of 149 spectra) as a func- 
tion of wavelength. The average and standard deviation of the 
total number of spectra are represented by a thin line. The 
derived template spectrum (/io, thick line) is identical to the 
average spectrum in the continuum, and should represent the 
undisturbed-wind profile in the variable Si IV line. In the lower 
panel the difference can be seen between the instrumental noise 
(estimated by ctq, thick line) and the observed variance (thin 
line) in the resonance lines due to wind variability. 



from DACs. The underlying assumption is that if the P Cygni 
profiles always return to the same flux level after passage of a 
DAC, a sufficient number of flux points per wavelength bin will 
be available (provided that the total number of spectra is large 
enough) to define the undisturbed-wind profile. 

We can compute /i and a 2 of these n highest flux values for 
each wavelength bin. We then assume that the "undisturbed" 
flux points belong to the wing of a normal distribution in any 
given part of the spectrum, which allows us to derive /io and <7q 
from /i and a 2 . In parts of the spectrum that do not vary intrin- 
sically (such as the continuum), the values obtained for /io (and 
a 2 ) should represent the average spectrum and the noise; this 
provides an independent check of the validity of our assump- 
tion. 



For our application we used the distribution function 4>(x) 
for a normal distribution given by: 



l 



<2it ctq 



I X fJ.Q \ ' 

\V2a J 



(1) 



The distribution of the n highest points is a truncated normal 
distribution, i.e. x running from w to oo, with w from Eq. (5). 
So we have for the average fi and variance a 2 of the n highest 
flux points: 



and 



,C x(j}{x)dx 
The ratio k is given by: 



with 



W - flQ 



(2) 



(3) 



(4) 



(5) 



y/2<ra ' 

After substitution of these expressions in Eq. (2) we can relate 
H to /i : 

2 



CT -z 

Ho = A* 1=- e 



'2-Kk 



(6) 



Thus, from jj, and ero we can compute /io; z is the inverse com- 
plementary error function of 2k (Eq. 4) and is computed itera- 
tively, using an adequate approximation for the complementary 
error function. Substituting /i from Eq. (6) in Eq. (3), we ob- 
tain a relation between (Tq and a 2 , the variance in the n highest 
points: 



2nk 2 



^2nk 2 + 2yfjrkze,- 



-2z 2 



(7) 



In principle, this solves the problem: Eq. (7) gives ao as a times 
a function of N and n only, and then Eq. (6) gives /i as a best 
estimate of the template. 

However, a is derived from a small number of points and 
has a large statistical uncertainty. Through Eq. (6) these uncer- 
tainties translate into noise in the template spectrum. Henrichs 
et al. (1994) have shown that for IUE spectra an empirical re- 
lation can be found between ctq and flux, which implies that 



we can derive ao from /io (see also Howarth & Smith 1995). 
This relation is based on much better statistics and gives more 
reliable results. Since /io is the quantity we want to derive, we 
have to estimate cto from /i, so that in principle the scheme has 
to be iterated. In practice, these iterations are barely needed. 
The method was tested numerically by selecting a given num- 
ber of highest flux points from a gaussian distribution (with 
known /xq and <tq) of randomly chosen points. The method was 
found to be extremely powerful: for k = 0.01 we predicted /io 
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Fig. 2. Left: the template spectrum (thick line) and a "representative" spectrum (thin line) including the Si IV resonance doublet 
are shown for the different stars. The units of the y-axis are in FN/s. Right: division of the observed spectrum by the template 
gives the quotient spectrum (thin line), clearly bringing out the DACs in the spectra. The best fit is illustrated as a thick line. 



within 2%. The higher the value of k, the better the obtained es- 
timate for it ( an d °oX but we are limited by the fact that only 
a few points per wavelength bin are not disturbed by intrinsic 
variability. 

In Fig. 1 we present, as an example, the Si IV resonance 
doublet of 68 Cyg. From the available 149 spectra we selected 
the 5 highest flux points per wavelength bin and computed the 
average fi and variance er 2 . We show //(A) and cr(A) by a dotted 
line in the upper and lower panel of Fig. 1 , respectively. The av- 
erage flux points and standard deviations of all 149 spectra are 



shown as a thin line. Note that the average of the 5 highest flux 
points fi(X) lies, as expected, systematically above the template 
/io(A), and that in the lower panel the scatter in a is large. The 
derived /zn (A) and <jq (A) are represented by a thick line. We see 
that the template /xo(A) and the average spectrum are identical 
in the continuum regions, but that the template deviates sig- 
nificantly from the average spectrum in the variable resonance 
line. In the lower panel of Fig. 1 the similarity between the es- 
timated cr (thick line) and the measured variance (thin line) is 
obvious in the wavelength regions outside the Si IV line. 
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Fig. 3. As Fig. 2: Template spectrum (thick line) and sample model fits (thick line) for the N V resonance doublet. 



Figs. 2 and 3 display the relevant parts of the derived tem- 
plate spectra. Shown are the P Cygni profiles including DACs. 
For comparison, a representative observed spectrum (thin line) 
is shown in the same panel. The quality of the template can 
be checked by searching for "emission" in the obtained quo- 
tient spectra (i.e. where the flux exceeds unity), which should 
be present if the template spectrum still contains additional ab- 
sorption with respect to some of the observed spectra. Another 
check is to compare the variance a 2 in the selected highest flux 
points in the variable line with the value of a 2 in a constant 
part of the spectrum with a comparable exposure level. For ex- 
ample, in the case of £ Per we were forced on these grounds to 
select only the two highest flux points. Although a logical con- 
sequence of our method is that a few points (e.g. at least one 
per wavelength bin for £ Per) in the quotient spectra will exceed 
unity, we found that these points are not randomly distributed 
among the quotient spectra, but are concentrated in only a few 
of them. This means that we can only improve on this by en- 
larging the sample of observed spectra. Fortunately, the "emis- 
sion" in the quotient spectra of £ Per is modest. For the other 
stars this problem turns out to be even less important. 

In order to investigate whether there exists such a thing as 
an undisturbed wind profile which represents the state to which 
the wind returns after the passage of a DAC, we constructed 
templates based on the spectra of the individual datasets and 
compared them. The most challenging case is that of £ Per, 
where the template we used is based on the two highest flux 
points out of a total of 124 spectra. As can be expected, the 
resulting templates are identical (within the noise) in the con- 
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Fig. 4. Comparison of templates in the Si IV region of £ Per 
constructed from 1994 data (thin line), 1991 data (dashed line) 
and all data prior to 1994 (thick line) used in the present analy- 



tinuum regions of the spectrum, but differ in the absorption 
troughs of the P Cygni lines. Of the templates based on the in- 
dividual datasets prior to 1994, the October 1991 template (36 
spectra) is closest to the template we used, but clearly contains 
more absorption in the blue doublet component (Fig. 4, Si IV 
doublet). We also constructed a template from a dataset of 70 
spectra of £ Per collected during a campaign lasting 10 days in 



October 1994 (Henrichs et al. 1998), which are not comprised 



in the sample studied in this paper. The October 1994 template 
resembles the template we constructed from all spectra until 
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1991 even more, but also here additional absorption is found at 
high velocities. This demonstrates that, given a sufficient times- 
pan covered by the dataset, our method is capable of finding 
the same underlying wind profile using different datasets; in 
the case of £ Per, however, only at low and intermediate ve- 
locities. For the latter star we have to conclude that in October 
1991 and October 1994 the high- velocity part of the wind has a 
larger optical depth compared to other campaigns, which could 
be due to long-term variability. 

2.2. Modeling DACs in quotient spectra 

The isolated DACs in a quotient spectrum, obtained after di- 
vision of an observed spectrum by the template, are modeled 
in the way described in Henrichs et al. (1983 see also Telt- 



ing & Kaper 1994). The DACs are assumed to be formed by 



plane-parallel slabs of material in the line of sight, giving rise 
to an absorption component with central optical depth r c and 
(Doppler) broadening parameter v t . The intensity of the com- 
ponent is described by: 



I(v) 



with the Gaussian profile function 

(v) = cxp - ( 



(8) 



(9) 



where v is the velocity with respect to the stellar rest frame and 
v c the Doppler displacement at line center. With the (reason- 
able) assumption that the displacement v c and broadening pa- 
rameter v t are identical for the absorption components in both 
lines of the resonance doublet, the DACs in both doublet com- 
ponents can be modeled simultaneously knowing the doublet 
separation v sp ij t and the ratio in oscillator strength. This leaves 
only three free parameters for each pair of DACs. These param- 
eters (v c , v t , and r c ) are determined by means of a x 2 -method 
in order to obtain the best fit. The signal-to-noise ratio in each 
point of the quotient spectrum is estimated using the empiri- 
cally determined dependence of this ratio on the flux (Paper I). 
The x 2 criterion then gives the f ormal errors in the derived pa- 
rameters (Telting & Kaper 1994). We fitted our model to spec- 



tral data within the velocity range —3000 to +2500 km s 
with respect to the principal (short-wavelength) doublet com- 
ponent. From these parameters we determine the column den- 
sity A^coi as sociated with the absorption components (Henrichs 
et al. 119831): 



m e c 



T c V t 



ne 2 f\ (1 + v c /c) 



(10) 



As an example, we show in Figs. 2 and 3 model fits (thick 
line) to quotient spectra which were obtained after division of a 
representative spectrum (thin line, left-hand panel) by the tem- 
plate (thick line). In most cases, a number of components were 
modeled simultaneously, with a maximum of 4 DACs in one 
doublet component. Furthermore, all fit parameters were kept 
free. As an initial condition we used the results of the previous 



f Per 07.5 III Si IV 
October 1991 

Window function 




0.5 1 1.5 2 2.5 3 3.5 4 

Frequency (cycles/day) 

Fig. 5. The window function obtained for the time series of 
£ Per in October 1991 (full line) is very clean thanks to the 
continuous coverage. For comparison, the cleaned power spec- 
trum, integrated over the red absorption component of the Si IV 
doublet, is shown as a dashed line. 



fit and modeled the spectra in a sequence determined by the 
time of observation. Sometimes we had to go back and forth 
through the dataset when it appeared that we had missed the 
first (weak) signs of the development of a new DAC at low 
velocity. In some cases we had to correct for bad normaliza- 
tion in the quotient spectrum by allowing for a very weak and 
wide (sometimes a few thousand km s -1 ) absorption compo- 
nent. This way we were able to determine the central velocity 
and width of the other components with reasonable accuracy. 
In these cases the derived column densities are, however, less 
accurate. 



2.3. Period-search analysis 

We also inspected the P Cygni lines for the occurrence of peri- 
odic variability. The spectral lines were normalized to the level 
of the surrounding continuum through division by a first-order 
polynomial. The period-finding algorithm we applied, consists 
of an ordinary Fourier transformation for non-equidistant tem- 
poral sampling, followed by a CLEAN stage in which the win- 
dow function, which is due to incomplete temporal sampling of 
the stellar signal, is iteratively removed from the Fourier spec- 
trum (Roberts et al. 1987 1. The continuous observations from 



space result in a very nice window function which does not 
include a peak at a period of one day (see Fig. 5 for an exam- 
ple). The window function was removed in 400 iterations with 
a gain of 0.2. The studied frequency range runs from 0.01 to 
10 cycle day -1 , with a frequency step of 0.01 c/d. A Fourier 
spectrum is computed for each wavelength bin (0.1 A). The 
power in the 2-dimensional periodograms resulting from the 
CLEAN analyses can, provided that the window function could 
be satisfactorily removed, be transformed back to amplitudes 
in continuum units (amplitude = 2 (power) 1 / 2 ), contrary to the 
summed power in the summed 1 -dimensional periodograms. 
More details on the used method and references can be found in 



Gies & Kullavanijaya (1988) and Telting & Schrijvers (1997) 
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It turns out that for several of our targets the derived period 
is close to, or sometimes even larger than the time span At 
covered by our observations. In these cases we can only con- 
clude that a long-term variation is present in the data, which 
may or may not be periodic. We have indicated these "periods" 
by putting them between square brackets (e.g. [6.2±1.8 d]). 

An outcome of the period-search analysis is the phase infor- 
mation. For a given period, the corresponding sinusoid's phase 
(in our figures in units of radians) is known in each wavelength 
bin. For example, periodic features moving from the red to the 
blue side of a spectral line will give rise to a shift in phase with 
wavelength. Owocki et al. (1995) used the observed phase di- 
agram of a periodic modulation in UV resonance lines of the 
BO. 5 lb star HD 64760 as supporting evidence for the occur- 
rence of corotating streams in the wind. In this particular ex- 
ample, however, the so-called "phase bowing" did not corre- 
spond to the evolution of the, though present, DACs, but to 
(sinusoidal) modulation s in fl ux at low and intermediate ve- 
locities (Fullerton et al. 1997 ). It would be very interesting to 
know whether these findings also apply to the here studied O- 
star wind variability. 



3. DAC behaviour in individual stars 



In this section we describe the results based on the modeling 
of DACs in the UV P Cygni lines and period-search analyses 
for each star separately. For some stars several time series are 
available, so that we can study the behaviour of DACs also on a 
longer (yearly) timescale. We have chosen to analyse in detail 
only those time series (from the sample presented in Paper I) 
that exhibit significant variability and cover a sufficiently long 
time interval. 

In most cases the Si IV resonance doublet was used to 
model the DACs; the other UV resonance lines (N V and C IV) 
are often saturated. The two main-sequence stars in our sam- 
ple, 15 Mon and 10 Lac, are variable in the N V (and C IV) 
profiles, but not in the Si IV doublet, which is too weak to show 
any wind features. For the supergiants 19 Cep and £ Ori we 
modeled the DAC behaviour in both the Si IV and the N V dou- 
blet. The subordinate N iv line at 1718 A of £ Per, HD 34656, 
and A Cep (for the latter also the He II line at 1640 A) varies 
in concert with the DACs in the Si IV doublet, but at low ve- 
locities only. In the wind profiles of a Cam we did not detect 
any significant variations, except for some marginal variability 
in the blue edge of the Si IV and C IV lines (Paper I). 

To facilitate a visual comparison, the figures showing the 
time dependence of the DAC parameters are plotted on scale. 
The scale of the time axis is identical in all figures. For a given 
star, the scale of the vertical axes is the same for each dataset. 

3.1. HD 24912 (£ Per) 07.5 III(n)((f)) 

We observed £ Per four times in the period 1987-1991. Fig. 6 
presents an overview of the DAC behaviour in the Si IV res- 
onance doublet and enables a visual comparison to be made 
between the different datasets. As discussed in Paper I, the vari- 



Fig. 9. £ Per October 1991: Two-dimensional periodogram of 
the Si IV doublet at 1400 A. The CLEANed Fourier spectrum 
obtained in each wavelength bin (horizontal axis) in the Si IV 
line is plotted as a function of frequency (in cycles day -1 , verti- 
cal axis). The power is represented in levels of grey (white cut: 
10~ 5 ; black cut: 10 -3 ). The mean Si IV profile is plotted in the 
bottom panel. For example, at a wavelength of 1388 A maxi- 
mum power is detected at a frequency of 0.5 c/d. An important 
result is that the Fourier spectrum is different at the highest ve- 
locities with a dominant peak at lower frequency. Also note that 
the power shows a gap at the asymptotic DAC velocity. 



ability pattern is characteristic for a given star, but shows de- 
tailed changes from year to year. The timescale of variability 
(i.e. the recurrence timescale of DACs) and the range in veloc- 
ity over which variability is observed, remain the same. The 
strength of the absorption components, however, varies from 
event to event. 

A minimum-absorption template was constructed from the 
124 Si IV spectra of £ Per, which was used to generate quotient 
spectra (Fig. 6). The isolated DACs were fitted to a model pro- 
file as described in section 2.2; in Figs. 7 and 8 the evolution 
of the DAC parameters (v c , Vt and r c ) are presented for the 
different sets of spectra. In the upper panels of Figs. 7 and 8 
the total equivalent width (integrated over velocity rather than 
wavelength) of the quotient spectra is given as a function of 
time. The error bars represent lcr-errors (cf. Telting & Kaper 



1994). We connected the points that we identify as belonging 



to the same DAC event. The criteria for selection of these points 
are based on the assumption of continuity in the three fit param- 
eters and the column density. 
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£ Per 07.5111 October 1991 




f Per 07.5111 1987-1991 

Si IV, C IV(x10), N IV(x20) 
September 1987 At=3.6d 



0.5 1 l.S 2 

Frequency [Cycles/day] 

Fig. 10. £ Per October 1991: One-dimensional power spectra 
(obtained by integrating a 2d-periodogram such as displayed 
in Fig. 9 over a given wavelength range) of the N V, Si IV, 
C IV, and N IV P Cygni lines. The periodic variations at low 
velocity are represented by a solid line, the variations at the 
blue edge by a dotted line. The 2-day period (0.5 c/d) appears 
in all lines at low velocity, while the blue edge seems to vary 
with approximately twice that period. 

The derived column density of a DAC first increases and 
then decreases; a maximum is reached at a velocity of about 
— 1500 km s _1 . The maximum in N co i is 5x 10 14 cm~ 2 . The 
other DAC parameters vary less strictly. The maximum width 
of a strong component is about 750 km s _1 just after its ap- 
pearance at low (~ —1000 km s _1 ) velocity. We have modeled 
narrow components at their terminal velocity having widths 
smaller than 50 km s _1 , i.e. nearly down to the spectral reso- 
lution of the IUE spectrograph, but not as narrow as the DACs 



observed with the GHRS onboard HST (Shore et al. |1993| ). The 
maximum central optical depth of a DAC that is observed for 
£ Per is r™ ax -1.2. 

Two different types of DACs can be discriminated in the 
Si IV doublet of £ Per: (1) strong absorption components that 
remain visible for more than two days (filled symbols) and 
(2) more short-lived, sometimes suddenly disappearing DACs 
(open symbols), preceding or following a strong component 
within half a day. 

For each dataset we performed a period analysis on the N V, 
Si IV, and C IV resonance lines and the N IV subordinate line. 
The Fourier analyses clearly detect periodicity, the cleanest pe- 
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Fig. 11. £ Per 1987-1991: A comparison of the (summed) Id 
power spectra calculated for the different datasets. The solid 
line represents the integrated power (x20) for the N IV line 
(—1000 to km s _1 ), the dotted line shows the power in the red 
doublet component of the Si IV profile (—1800 to km s _1 ), 
and the dashed line gives the power (x 10) in the blue edge of 
the C iv doublet (-2950 to -2650 km s" 1 ). 



riodogram being obtained for the October 1991 dataset (see 
Fig. 9 for the Si IV results), when £ Per was covered for the 
longest period of time (At = 4.4 d; consequently, all periods 
derived from this dataset that are close to or longer than 4.4 day 
are not well determined and will be printed between square 
brackets). In each wavelength bin a CLEANed Fourier spec- 
trum is calculated and plotted as a function of frequency (in 
cycles day -1 ), the power being represented in levels of grey. 
Periodic variability is detected in the blue-shifted absorption 
troughs of the Si IV doublet (also in the other studied lines). 
The dominant period is found at 2.0±0.2 day, almost over the 
full width of the blue-shifted absorption profile. As a conser- 
vative estimate of the accuracy of the period we take the half- 
width of the peak in the (summed) power spectrum. The 2-day 
period is clearly reflected by the change in total EW with time 
(Figs. 7 and 8). 

At higher velocities (>1800 km s _1 ), however, a peak is 
found at a longer period: in the Si IV line at [5. Oil. 5] day, 
which is longer than the length of the observing run and there- 
fore unreliable. In the other lines (N V and C IV) also a longer 
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Fig. 12. £ Per 1987 and 1991: Phase (in radians) of the 2.0- 
day period as a function of position in the line. The phase is 
only plotted when the summed power is larger than 0.0001. 
The 1987 points are displayed as crosses, the 1991 points as 
filled circles. The two phase diagrams are matched at an arbi- 
trarily chosen velocity for comparison. Phase bowing is clearly 
observed in the Si IV doublet. 



period is detected at the blue edge (between —2950 and —2350 
km s" 1 ) of the profiles, at [3.8±1.0] and [3.6±0.9] day, respec- 
tively (Fig. 10, dotted lines). These (saturated!) resonance lines 
also include the 2-day period at lower velocities (ranging from 
about —2000 to km s _1 with respect to the rest wavelength 
of the red doublet component). The longer period is in all three 
cases consistent with a period of 4 days, i.e. twice the period 
observed at lower velocities. The subordinate N IV line only 
exhibits the 2.0±0.3 day period, at velocities between —1000 
and km s _1 . 

The 2-day period is also present in spectra from previous 
campaigns. In Fig. 1 1 we compare ld-power spectra from dif- 
ferent lines for the period 1987-1991. The solid line shows the 
power spectrum of the N IV line, integrated over the absorp- 
tion part of the line from —1000 to km s _1 . The campaigns 
covering the shortest period of time (1988 and 1989) result in a 
higher weight to a 1-day period (1.0±0.1 d). The evidence for 
a longer period in the blue edge of the profiles, such as found in 
1991, is not very strong. This is not surprising given the short 
time span covered by these observations. 

Why do the observed variations at higher velocities exhibit 
a longer period, about twice as long as at lower velocities? If 
our interpretation is correct, this factor of two difference in pe- 



riodicity is due to the peculiar DAC behaviour in the wind lines 
of £ Per. In Paper I we pointed out that in this star strong DACs 
reach different asymptotic velocities, leading to crossings of 
successive components. We found that the DAC's asymptotic 
velocity alternates between 2050 km s _1 and 2275 km s _1 
(in Figs. 7 and 8 the high-velocity DAC is indicated by filled 
squares). Therefore, only one out of two DACs would reach 
the far blue edge of the profile resulting in a frequency of once 
every four days at the highest velocities. 

It turns out that if one wants to understand the edge vari- 
ability, a detailed knowledge of the DAC behaviour is essential. 
The recognition that the asymptotic velocity of a DAC can have 
two different values has important consequences for the use of 
this diagnostic to determine the terminal velocity of a stellar 
wind (Henrichs et al. 1988, Prinja et al. 199C ). We assume here 
that the highest value corresponds to Voo, Furthermore, the al- 
ternating behaviour of DACs in the wind of £ Per suggests that 
a full cycle includes two strong DAC events, i.e. 4 days. 

Fig. 12 shows the phase of the 2.0-day period as computed 
by the Fourier analysis. Only the results are plotted for the 1987 
(crosses) and 1991 (filled circles) datasets; the 1988 and 1989 
datasets give similar results, although there a 1-day period is 
dominant (the phase diagrams are, however, almost identical 
to those of the 2-day period). The phase (in radians) is plotted 
as a function of position in the line (measured in velocity with 
respect to the blue doublet component) for the three resonance 
lines and the N IV subordinate line. The phase has only a phys- 
ical meaning if at the given velocity the period is detected. We 
plot the phase for points with power exceeding 10~ 4 , which is 
just above the level of the continuum in the power spectrum 
(Fig. 10). The error in phase does not follow straightforwardly 
from the Fourier analysis; a possibility is to bin the data and 
use the spread in phase as an indication of the error. We did not 
include error bars in the phase diagrams, but the spread in the 
points indicates that the error on the phase must be small. 

The Si IV lines clear ly dem onstrate the s o-calle d phase 
bowing (cf. Owocki et al. 1995 Fullerton et al. 1997 ), indica- 
tive of curved wind structures (like corotating interaction re- 
gions) moving through the line of sight. A maximum in phase is 
reached at about —1200 km s _1 in the Si IV profile. In the CIR 
model, this velocity corresponds to the velocity of the material 
in the spiral-shape CIR first moving out of the line of sight. The 
maximum phase difference A<fi (measured between w —1000 
and —2000 km s _1 ) is about one radian, which is comparable 
to that observed for the B sup ergiant HD 64760 (A</> w 0.6 
7T radians, Fullerton et al. 1997). Note that a decrease in phase 
towards more (less) negative velocities corresponds to an ac- 
celerating (decelerating) feature moving upwards in time (as in 
the grey-scale plot of e.g. Fig. 6). Both datasets show a similar 
picture; for comparison, an arbitrary constant was added. The 
N IV line shows the rising phase only at low velocities, consis- 
tent with the behaviour of the Si IV (and N V) doublet. Satu- 
ration of the C IV and N V lines causes the absence of points 
bluewards of —1000 km s _1 . At low velocities the 1987 points 
seem to indicate a faster increase in phase. Another small dif- 
ference between 1987 and 1991 is seen at high velocities and in 
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HD34656 07 11(f) February 1991 



X Ori Si IV DAC fits November 1992 
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Fig. 14. HD 34656 February 1991: the top three panels present 
ld-periodograms illustrating the periodic variability in the N V, 
Si IV, C IV, and N IV P Cygni lines. The solid lines show the 
power summed over the low-velocity (<1500 km s _1 ) part of 
the profiles, the dotted lines represent the high-velocity part. A 
[4.5-day] period is only found at high velocities, a 1- and/or 2- 
day period dominates at low velocity. The bottom panel shows 
the phase (in radians) of the [4.5] (filled circles), 1.8 (open 
squares), and 1.1 -day period (crosses) as measured from the 
Si IV line (power >1CT 4 ), only for the blue doublet compo- 
nent (An = 1393.755 A). 



the blue edge. In 1987 the edge variability is very pronounced 
and leads to a signal at the 2-day period around —2500 km s _1 
in the resonance lines as mentioned earlier. 



3.2. HD 34656 0711(f) 

We observed the 07 supergiant HD 34656 only once; the 27 
spectra obtained during 5.1 days in February 1991 show the 
rapid recurrence of several DACs appearing at low velocity 
(about —500 km s _1 ). In Fig. 13 the fit parameters are shown 
for the DACs in the Si IV line. The DAC behaviour in this star is 
quite complicated: DACs appear about once every day (observ- 
able as regular humps in EW), but around day 10 the strength 
of the absorption at —1500 km s _1 increases rapidly. It appears 
as if an additional (variable) absorption component is present 
at —1500 km s _1 at all times. This makes it difficult to assign 
points to specific DAC events. 




Time (HJD 



7 

2448930) 



Fig. 15. As Fig. 7: the parameters of the migrating DAC in the 
Si IV profile of the 08 HI((f)) star A Ori in November 1992. 



The column density of this additional component suggests 
that it belongs to a separate DAC event, reminiscent of a per- 
sistent component as observed in the wind lines of A Ori 
and 15 Mon (see below). At day 9.8 the 7V co i of this compo- 
nent reaches a minimum of 4xl0 13 cm~ 2 and then starts to 
strengthen rapidly to a maximum of 2.5 x 10 14 cm~ 2 . We sus- 
pect that at day 9.8 a new "DAC" replaces the former one at 
— 1500 km s _1 and remains at that velocity until day 12. This 
would mean that the absorption component at —1500 km s _1 
follows a different cycle than the "rapid" DACs appearing once 
every day. The variations might be interpreted as caused by 
two DAC patterns superimposed on each other: a "slow" pat- 
tern consisting of DACs with an asymptotic velocity of — 1500 
km s _1 , and a "fast" pattern of DACs that reach higher veloci- 
ties (<~ —2150km s _1 ). The rapid DAC pattern is also (weakly) 
present in the N IV line, the DACs can be traced down to a ve- 
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locity of —200 km s _1 . The origin of the slow pattern is not 
clear; it might reflect a variation in wind structure on a larger 
spatial scale than the structures causing the DACs belonging to 
the rapid pattern. 

Period analyses of the P Cygni lines of HD 34656 also 
indicate a bimodality in the DAC behaviour. Fig. 14 presents 
ld-periodograms obtained after summing the 2d-periodograms 
over a given wavelength interval. The solid lines show the 
power detected at low velocity (<1500 km s -1 ), the dotted 
lines the periodic behaviour at the blue edge. The dominant pe- 
riod found in the N IV line is 1 . 1 ±0. 1 d and corresponds to the 
DAC recurrence timescale in the "rapid" pattern. This period 
is present in the low-velocity part of all four lines. A period of 
1.8±0.3 d is detected in all lines (at low velocity), and might 
be related (alias) to the 1-day period variation. At high velocity 
(dotted lines) a longer period dominates: [4.5±1.8 d] (in the 
case of Si IV). This period might relate to the slow pattern (see 
also the corresponding changes in EW, Fig. 13). 

The phase (in radians) of the 3 detected periods is plotted 
as a function of position in the Si IV line in the bottom panel 
of Fig. 14. The phase is plotted for the blue doublet component 
only. In this line, most power is measured for the 4.5-day period 
at velocities bluewards of —1500 km s _1 . For all periods, the 
phase shows a declining trend towards higher velocities, but no 
obvious phase bowing, such as detected for £ Per, seems to be 
present. 

3.3. HD 36861 (X Ori) 08 III((f)) 

We observed the 08 III((f)) star A Ori during the November 
1992 campaign and obtained 27 spectra (At = 5.0 d). The 
Si IV resonance lines of A Ori consist of a sharp photospheric 
absorption component, slightly asymmetric to the blue due to 
wind contamination, plus an absorption component at —2000 
km s _1 (Paper I). Also the N V and C IV resonance doublets 
exhibit this persistent absorption feature at —2000 km s -1 , in 
both doublet components. We modeled the quotient Si IV spec- 
tra of A Ori (Fig. 15) and found one DAC accelerating from 
— 1700 to —2000 km s _1 , i.e. reaching the position of the per- 
sistent component. The occurrence of the DAC is reflected by 
the rise in Si IV EW around day 5.7, shown in the upper panel 
of Fig. 15. We covered one DAC event, resulting in a lower 
limit of about 5 days to the DAC recurrence timescale. These 
observations suggest that the persistent component is regularly 
"filled up" by new absorption components when they reach the 
terminal velocity of the wind (which can, if this interpretation 
is valid, be determined from the position of the persistent com- 
ponent). 

Compared to DACs in other stars in our sample, the mea- 
sured column density is low, about 5 x 10 13 cm -2 in Si IV. The 
central optical depth does not exceed 0.3 while the width of the 
component is smaller than 250 km s _1 . Both in the N V and 
in the Si IV line a period around 4 days is detected (between 
1500 and 2000 km s^ 1 ): [3.4±0.8 d] and [4.7±1.4 d], respec- 
tively. Since only one DAC event was observed, the physical 



significance of these periods is doubtful. It might just reflect 
the presence (and subsequent absence) of a DAC. 

3.4. HD 37742 (( Ori) 09.7 lb 

Moving DACs are found in both the N V and the Si IV reso- 
nance doublet of the 09.7 lb star ( Ori. The results from our 
fit procedure are shown in Fig. 16. For this star it is very dif- 
ficult to identify the different DAC events. Our reconstruction 
is as follows: At day 5.4 and 7.1 the appearance of a DAC is 
registered in both lines, with the difference that the first com- 
ponent accelerates faster than the second one, and reaches a 
velocity of —2000 km s _1 ; the second component ends at a 
much lower velocity (—1250 km s _1 ), similar to the compo- 
nent present since the beginning of the observations. At day 
8.7 a new component develops, consistent with a recurrence 
timescale of about 1.6 days. (Alternatively, bearing in mind the 
v c versus time plot (Fig. 16), the points at low (—500 km s _1 ) 
velocity at day 5 could be connected to the event starting at day 
7 (open squares); the filled squares would then be a continua- 
tion of the open and filled circles, respectively.) 

The Fourier analysis produces a peak in power at a 
1.6±0.2 d period; most of this power is concentrated in the 
velocity range —1500 to —2150 km s _1 . The highest peak in 
the power spectrum is found at a period of [—6] days, which is 
longer than the 5.1 d span of the time series. Ha observations 



of this star (Kaper et al. 1998) indicate a 6-day period as well, 



so that this period might be real. The 1.6-day period is about a 
quarter of the 6-day period. 

The properties of the DACs in the two different lines can 
be compared. The central velocities of the modeled compo- 
nents are similar in both lines. The other DAC parameters 
show the same trend in both the N V and the Si IV line. The 
column densities differ by a factor of about two: maximum 
A col -1.1 xlO 14 and 0.6xl0 14 cm" 2 for N V and Si IV, re- 
spectively. The similarity of the DAC behaviour in the N V and 
Si IV lines supports the common view that the variable DACs 
reflect changes in wind density (and/or the velocity gradient), 
rather than in the ionization structure of the stellar wind. 

3.5. HD 47839 (15 Mon) 07 V((f)) 

Two migrating DACs are found in the N V doublet of the 
07 V((f)) star 15 Mon. The first DAC is visible from the start 
of the campaign, the second one appears just before day 9 (Fig. 
17) and is very narrow (w f < 150 km s _1 ). Such a narrow width 
is also observed for DACs in P Cygni lines of A Ori and 10 Lac. 
Because we have only 20 spectra in total, it is difficult to con- 
struct a good template; this explains why the EW of fits to some 
quotient spectra is negative (not shown). Furthermore, the vari- 
ations in the N V profile have a small amplitude, which makes 
the quality of the template an even more important factor. The 
first DAC accelerates up to —2250 km s — 1 , which is consid- 
erably higher than the central velocity (—1900 km s _1 ) of the 
persistent component in the N V (and C IV) profile. The sec- 
ond component, which appears at a velocity of — 1700 km s , 
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seems to join the persistent component. The column density of 
the DACs is less than 10 14 cm -2 , but this is uncertain due to 
the poor normalization. 

On the basis of this dataset, only a lower limit of 4.5 days 
can be set to the DAC recurrence timescale. The Fourier analy- 
sis indicates a [6-day] period (power integrated from —1500 to 
—2350 km s -1 ), but the length of this period exceeds the 5.3 d 
observing period. 

3.6. HD 203064 (68 Cyg) 07.5 III:n((f)) 

The August 1986 dataset of the 07.5 I II:n((f )) star 68 Cyg 
has been analyzed by Prinja & Howarth ( 1988 ). In Fig. 18 we 
show the results of profile fits to this set of 33 spectra. Our re- 
sults are compatible with those of Prinja & Howarth. In Fig. 19 
the DAC parameters derived for the September 1987 and Oc- 
tober 1988 campaigns are presented. The central velocities of 
DACs in the September 1987 spectra were presented by Fuller- 



68 Cyg 07.5111 1986-1991 



ton et al. (|1991j). The DAC fits of the October 1989 and October 
1991 datasets are shown in Fig. 20, the latter dataset having the 
longest time coverage with good time resolution (5 days); we 
describe these observations first. 

The Fourier analysis results in a pronounced peak in the 
power spectrum at a frequency of 0.73±0.19 d~ x , correspond- 
ing to a period of 1.4±0.2 day (Fig. 21). Assuming that this 
period is the DAC recurrence time scale, we can identify two 
different series of DACs moving through the Si IV profile (Fig. 
20): series A, which appears at JD 52.8, 53.9, 55.5, 56.5, and 
series B, appearing at JD 53.4, 54.7, 56.0 (the sampling time 
of the UV data is 0.1 day). The time lag between the two se- 
ries is about half a day, series A followed by B. Furthermore, it 
turns out that DACs belonging to different series, reach differ- 
ent asymptotic velocities: series A accelerates to a velocity of 
about —2450 km s _1 , which is ~250 km s _1 higher than the 
velocity reached by series B. Thus, like in the case of £ Per, the 
DACs asymptotic velocity shows an alternating behaviour. An 
important difference is, however, that in £ Per the time span be- 
tween two successive DACs with different asymptotic velocity 
(i.e. belonging to different series) is equal to the DAC recur- 
rence timescale, while in 68 Cyg the time lag between series A 
and B is about 0.5 days, i.e. not equal to the DAC recurrence 
timescale (1.4 days). 

Simultaneous Ha observations indicate that incipient Ha 
emission is observed prior to the appearance of DACs belong- 



ing to series A (Kaper et al. 1997), and not in between series A 



and B. Alternatively, each DAC event in 68 Cyg might consist 
of two separate components, the event starting with an absorp- 
tion component from series A, followed half a day later by one 
from series B. 

When comparing the different datasets, typical DAC pat- 
terns are recognized (see also the figures in Paper I). In Paper I 
we pointed out the remarkable similarity between the DAC pat- 
tern in the Si IV time series obtained in 1986 and 1988. In the 
figures we labeled the DACs showing similar behaviour with 
the same symbol. Like for £ Per, the DAC behaviour suggests 
that the full cycle of variability is twice the DAC recurrence 
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Fig.21. 68 Cyg 1986-1991: A comparison of the (summed) 
Id power spectra calculated for the different datasets. The 
solid line represents the summed power for the blue absorp- 
tion component of the Si IV line (—2100 to —600 km s _1 ), the 
dashed-dotted line that for the red doublet component (—2100 
to —150 km s _1 ). As expected, both lines follow the same 
trend. The long-dashed line demonstrates the periodicity found 
at the high-velocity edge of the Si IV profile (-2750 to -2100 
km s _1 ). The dotted line gives the integrated power spectrum 
obtained for the red doublet component of the N V profile 
(—1875 to —675 km s _1 ). The short-dashed line shows the 
power spectrum in the blue edge of the C IV doublet (—2800 to 
-2550 km s" 1 ). 



timescale, which is 2.8 days. In exceptional cases (as, for ex- 
ample, in the time series of September 1987 demonstrated by 
the absence of filled squares) a DAC seems to be missing. This 
is another indication that the general DAC behaviour is very 
regular over the years, but that small though significant differ- 
ences do occur. We note that with our fit procedure it is difficult 
to resolve two overlapping absorption components of compa- 
rable strength, and that such double features will sometimes fit 
as one single component, which causes unavoidable disconti- 
nuities in the time sequence of the fit parameters. 

The total DAC EW (upper panels) varies roughly with the 
1.4 day period. The column density of the DACs reaches a 
maximum when the component is half-way its acceleration 
(i.e. at a central velocity around —1750 km s -1 ), and then de- 
creases with time. The strongest DACs reach a column density 
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Fig. 22. As Fig. 12: Phase diagrams for the 1.4-day period (0.73 
c/d) detected in the Si IV lines of 68 Cyg (phase in radians). The 
datasets from different years give very similar results. Phase 
bowing is not convincingly detected, this might be due to the 
absence of variability at low velocities. 



of 3xl0 14 cm" 2 . Some components remain visible for two 
days. The maximum central optical depth t c is about 1, and 
the width of some DACs exceeds 1000 km s _1 when they first 
appear at low velocity. 

Fig. 21 presents the (summed) Id power spectra obtained 
for different UV P Cygni lines in the period 1986-1991. The 
power diagrams are quite complicated: at a given velocity sev- 
eral periods are detected, which gradually change from one ve- 
locity bin to the other. Therefore, excluding the October 1991 
dataset, the summed Id power spectra contain several rela- 
tively broad peaks. It is not possible to transform the summed 
power back to a variability amplitude in continuum units. In all 
years, a period of 1 .4 days is present, the campaigns covering 
the shortest period of time (1988, 3.1 d, and 1989, 2.6 d) pro- 
vide less significant results. The October 1991 dataset includes 
a 2.4±0.4 d period in the high-velocity edges of the profiles, 
but it remains uncertain whether this period could be identi- 
fied with the suggested 2.8 d full cycle of wind variability. The 
other datasets produce power at periods around 3 days, but the 
significance of any peak in power at a period close to 2.8 days 
is not high. 



The phase diagrams for the 1.4-day period (phase in ra- 
dians) are very similar from year to year (Fig. 22). The dif- 
ferences are mainly due to changes in the amplitude of the 
variability (the phases are only plotted for velocity bins with 
power >10~ 4 ). The September 1987 dataset shows the phase 
behaviour over the largest velocity domain. Phase bowing is not 
convincingly detected, but this might be due to the absence of 
variability at lower velocities (also the N IV line is not variable, 
contrary to the N IV line in £ Per). 

3.7. HD 209975 (19 Cep) 09.5 lb 

The 09.5 lb star 19 Cep provides a very clear picture of 
the DAC phenomenon. Because the characteristic timescale of 
variability is much longer (about 5 days) than in some other 
well-studied stars, the behaviour of DACs is easily recognized. 
In Fig. 23 we present the time evolution of the different fit pa- 
rameters for the Si IV and N V spectra obtained in August 1986. 
The model fits to quotient spectra of both resonance lines give 
comparable results. In August 1986 we detected three DAC 
events. A very strong absorption component appears at low ve- 
locity (v c = —300 km s _1 and vt = 1000 km s _1 ) at day 

6.8. The column density of this component reaches a maximum 
of 6.8x 10 14 cm~ 2 at a central velocity of -1450 km s _1 (in 
Si IV). This is the strongest DAC we encountered in our col- 
lection of O-star spectra. The maximum central optical depth 
exceeds 2 at the peak column density. The asymptotic velocity 
of the two weak components present during the first half of the 
campaign is —1900 and —1750 km s -1 for the first and second 
component, respectively. 

The time series obtained in September 1987 (At = 3.2 d) 
and October 1988 (At = 3.8 d) spanned only a few days and 
did not reveal the clear evolution of a DAC (Paper I), which can 
be expected on statistical grounds when the DAC recurrence 
timescale is on the order of five days. It could also indicate that 
the strength of the DACs varies from year to year. In October 
1991 (Fig. 24) a weak component is present at its final veloc- 
ity of —2050 km s _1 ; at day 3.2 a new DAC appears at low 
velocity. The maximum column density reached by this com- 
ponent is 1.9xl0 14 cm -2 and r™ ax = 0.7. At day 5.9 a new 
DAC occurs with rapidly growing strength. The time interval 
of 2.7 days is short if we realize that in August 1986 the DAC 
recurrence timescale would be estimated to be on the order of 
5 days. The November 1992 observations (Fig. 25) again show 
the development of a DAC, a fairly strong one with maximum 
N co i of 3.5 x 10 14 cm -2 . At the end of this run (day 9.5) a new 
DAC appears in the Si IV line, from which we again would con- 
clude that the recurrence timescale is about five days. The 2.7 
days observed in October 1991 might correspond (within the 
errors) to half this period. 

The Fourier analysis of the October 1991 campaign does 
not produce a significant peak in power at a period close to 
2.5 days. All three studied datasets include a long period, al- 
ways close to or even longer than the length of the campaign: 
[7.0±2.3 d] (1986, At = 6.7 d), [6.3±2.1 d] (1991, At = 
4.2 d), and [4.5±1.9 d] (1992, At = 5.6 d) in the Si iv profile. 
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The latter period should get the highest weight. Fourier analysis 
of the N V line produces similar results. The obtained periods 
are consistent with the period of 5 days suggested above. Also 
the Ha observations that were carried out simultaneously with 
the October 1991 IUE campaign indicate a 5-day period (Kaper 
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et al. 1997). 
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Fig. 28. A Cep 06 I(n)fp 1989-1991: Summed Id power spec- 
tra showing the periodical variability in the wind lines. The pe- 
riodicity in the Si IV line is measured in three intervals: the red 
and blue doublet component (solid line: —600 to —1600 and 
dotted line: —200 to — 1400 km s _1 , respectively) and the high- 
velocity edge (long-dashed line: —2200 to —2650km s _1 ). The 
C IV edge is represented by the short-dashed line (—2200 to 
—2750 km s _1 ). The periodograms obtained for the subordi- 
nate N IV (long-dashed-dotted line: -450 to -1500 km s _1 ) 
and He II (short-dashed-dotted line: -200 to -350 km s _1 ) 
lines are shown as well. 



Although the 06 I(n)fp star A Cep was monitored during 
six campaigns, we present here the DAC modeling results only 
for the last campaign in October 1991. The saturation of the 
UV resonance lines makes the modeling of DACs very difficult. 
We performed a Fourier analysis on the October 1989, Febru- 
ary 1991, and October 1991 datasets. For the October 1991 ob- 
servations of the Si IV doublet we succeeded in measuring the 
central velocity of DACs and were able to determine the recur- 
rence timescale. The velocity domain from —1600 to —2000 
km s _1 is, however, difficult to model, since the quotient spec- 
tra fluctuate with large amplitude due to the division by small 
numbers and the low signal-to-noise. The DACs in the Si IV 



Oct 1991 
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Fig. 29. Phase diagrams obtained for the 1 .4-day period (filled 
circles) in the Si IV doublet of A Cep; only the red doublet 
component is shown (Aq = 1402.77 A). 



doublet (Fig. 26) behave similarly as observed in other O stars: 
the column density increases after the appearance of a DAC and 
reaches a maximum (in this case 2.5 x 10 14 cm -2 ) when a cen- 
tral velocity of about —1500 km s _1 is reached, although this 
is rather difficult to determine because of the problems men- 
tioned earlier. In every spectrum, an absorption component is 
fitted at a velocity close to —2100 km s _1 . Probably, this com- 
ponent is at the terminal velocity of the wind and the new DACs 
merge with this component when they approach their asymp- 
totic velocity (i.e. Voo). Due to the above mentioned limitations 
we cannot demonstrate this explicitly. 

Signs of DACs are found in quotient spectra of the subordi- 
nate He II (1640 A) and N IV (1718 A) lines (Fig. 27). The first 
traces of these DACs occur at very low velocity: —200 km s _1 . 
A clear difference with the behaviour of the N IV line in £ Per is 
that the variations in the He II and N IV lines of A Cep show the 
rapid acceleration of DACs such as observed in, e.g., the Si IV 
resonance doublet. The variations due to DACs in the N IV line 



of the 04 supergiant £ Pup are very similar (Prinja et al. 1992 ). 
We conclude that also in A Cep the observations indicate that 
DACs are formed already close to the star. Simultaneous Ha 
observations support this conclusion (Kaper et al. 1997 ). 

The results of the Fourier analysis are presented in Fig. 28. 
The summed Id power spectra obtained for the October 1989 
and February 1991 campaigns are shown for comparison. The 
October 1991 dataset shows a period of 1.4±0.2 d (the solid 
line gives the power in the red absorption component of Si IV 
between —600 and — 1600 km s -1 ); this period is also observed 
in February 1991, but then the dominant period is [4.3±1.2 d]. 



L. Kaper et al.: Long- and short-term variability in O-star winds 



15 



Visual comparison of the time series (Paper I) suggests that 
the variability timescale in the February 1991 data is indeed 
longer than in October 1991 (but the sampling of the dataset is 
worse). The N IV profile (long-dashed-dotted lines) produces 
the 1 .4 d period ( 1 .3±0. 1 d) too, the He II line vaguely indicates 
a possible harmonic at 0.72±0.03 d. The October 1991 data 
also show evidence for a period close to [5] days (although 
this is longer than the 4.5 d coverage). The longer periods are 
predominantly found in the high-velocity edges of the Si IV 
(long-dashed lines, -2200 to -2650 km s _1 ) and the C IV 
(short-dashed lines, —2200 to — 2750 km s -1 ) PCygni profiles. 

Although the October 1989 campaign lasted only 2.6 days, 
it appears that the signature of the variability is different com- 
pared to October 1991 but rather similar to that in the February 
1991 dataset. The dominant period is at [2.5±0.8 d] in Si IV 



and [2.2±0.7 d] in the C IV edge. Henrichs (1991 ) showed that 



the optical He II 4686 A line exhibits (wind) variations that are 
in phase with the modulation of the C IV edge, providing com- 
plementary support for the significance of this 2-day period. 

How can we reconcile the observed periods with each 
other? If 1.4 d were the DAC recurrence time, then 2.8 d is 
two-times and 4.2 d three-times this period. These periods are 
all within the error bars of the observed ones. We will come 
back to this point in section 4.6. 

The phase diagrams (Fig. 29) are difficult to interpret. The 
phase (in radians, for power >5x 10~ 5 ) related to the 1.4-day 
period in October 1991 shows a decline in the red doublet com- 
ponent of the Si IV doublet consistent with the DAC accelera- 
tion. In February 1991 the phase behaviour is not clear and 
in 1989 the 1.4-day period is not detected (except in the blue 
edge). There is no clear evidence for phase bowing in A Cep. 

3.9. HD 214680 (10 Lac) 09 V 

We modeled the DACs in quotient spectra of the 09 V star 
10 Lac which was observed in November 1992 (At = 5.0 d). 
DACs are present in both the N V and the C IV resonance dou- 
blet; here we show the results from the profile fitting of the 
N V doublet only, because the DACs are most prominent in 
this line. The strongest component (Fig. 30) appears at day 7.3 
and is preceded by a DAC starting 1 .4 day earlier, which is not 
detected in the C IV doublet. The maximum central velocity 
observed for both components is — 990 kms -1 . Maximum col- 
umn density of the strongest component (A™* = 10 14 cm -2 ) 
is reached when v c is about —875 km s -1 . The time series 
is too short to determine the recurrence timescale with confi- 
dence. The characteristic timescale of variability is, however, 
longer than 4 days if the modeled absorption components be- 
long to the same event, as was observed in the case of £ Per 
and 68 Cyg. Another argument in favour of a timescale signif- 
icantly longer than 1.4 day is that the acceleration timescale, 
in this case at least several days, is usually comparable to the 
DAC recurrence timescale. 

The Fourier analysis applied to the N V time series of 
10 Lac results in a peak at [6.8±2.3 d] (and not at 1.4 d), but 
this period is longer than the time coverage of the data. A more 



recent series of IUE spectra, covering about 20 days, was ob- 
tained in August 1995 (Henrichs et al., in prep.). In this se- 
ries, the N V and C IV doublet clearly show two complete DAC 
events, with a recurrence timescale of about 7 days. A Fourier 
analysis yields a period of 6.8±1.0 d, consistent with the one 
detected in the much shorter November 1992 dataset. 

4. Comparing DAC behaviour in different stars 

Some of the general characteristics of the DAC phenomenon, 
such as the narrowing of DACs when they accelerate through 
the profiles, their always increasing speed until they reach an 
asymptotic velocity, and their cyclical appearance are known 
for quite some time. In the following, we will evaluate the DAC 
properties in O-star wind lines as encountered in our observa- 
tions. 

4.1. The DAC parameters 

For the nine stars in our sample that show the presence and evo- 
lution of DACs in one or more of the time series, we have quan- 
titatively modeled the DAC behaviour. We have demonstrated 
that the shape of the individual absorption components can be 
very well described by an exponential Gaussian. Table 2 sum- 
marizes the main results from the fitting procedures. For each 
time series the extreme values of the fit parameters are listed. 
Obviously, when more than one DAC event was covered in the 
time series, the extremes do not necessarily correspond to the 
same event. 

The central velocity asymptotically reached by a DAC at 
the end of its track is a measure of the terminal velocity Voc of 
the stellar wind (Tab. 1). This velocity is typically 10-20% less 
than the maximum velocity in the steep blue edge of saturated 
P Cygni lines, but about equal to the velocity (fbiack) corre- 
sponding to their most blue-shifted point of saturation (Prinja 
et al. 199C). In radiation-driven winds is directly related 



to the stellar escape velocity which does not change on short 
timescales. In stars for which we collected more than one time 
series, the measured maximum values of Vc Symp are the same 
within 10%. We note, however, that in some cases we found 
systematic differences in vl symp from event to event (see sec- 
tion 4.2). The data do not suggest a relation between Wc' symp of 
a given DAC event and its (maximum) strength (r™ 2 " and/or 

The minimum central velocity, v™ ln , at which DACs can 
be detected is typically in the range 0.2-0.4 v^. However, 
for the main-sequence stars in our sample (A Ori, 15 Mon, 
10 Lac) vf [n is of the order 0.7-0.9 Also ( Oph (Q 9.5 V , 
vf n -1100 km s~\ Voo= 1480kms^ 1 ,Howarthet al. |l993^ 
with 0.74 Voc and ( Pup (04I(n)f, vf m -600 kmr 1 ,^ 
2520 km s~\ Prinja et al. |1992[) with 0.24 Voo fit this picture. 



There is no obvious link between the projected rotational ve- 
locity v sin i and the lowest velocity at which DACs are ob- 
served. Given the large initial width of a DAC (sometimes up to 
0.5 Voo), though usually small initial optical depth, wind vari- 
ability can be traced down to zero velocity in several cases (see 
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also Paper I). The difference in u™ n between main-sequence 
and evolved O stars might be due to the relatively weak winds 
of the main-sequence stars. Also the Ha line of main-sequence 
stars does in general not show any signs of wind variability 



(Kaper et al. 1998) 



The maximum optical depth (r™ ax ) encountered for these 
absorption features is around one. The strongest DACs are 
found in the supergiant and giant spectra, although £ Ori has 
relatively weak absorption components for a supergiant and 
£ Per relatively strong components for a giant. The DACs in 
main-sequence stars have the smallest optical depths. The min- 
imum DAC widths (v™ m ) we measured are several tens to a 
hundred km s _1 . The resolution of the SWP camera aboard 
IUE is about 25 km s _1 . 

The DAC maximum column density varies from event to 
event, but in general -ZV™ X (i.e. the highest value of jV co i in 
a given time series, often containing several DAC events) is 
largest for the supergiants and smallest for the main-sequence 
stars in our sample. The velocity at which a maximum in col- 
umn density is reached is always in the range 0.7-0.9 Voo, 
i.e. about half-way the evolution of a DAC (Figs. 31 and 
32). For ( Ori and 19 Cep the DACs were fitted in both the 
Si IV and the N V doublet. The results are consistent; the 
measured N co \ is about a factor 2 higher in the N V doublet. 
With reference to Eq. 10, the difference in oscillator strength 
(/(Si IV) //(N v) = 3.4) is the main factor. 
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Fig. 32. The variation of DAC column density as a function of 
central velocity in the Si IV resonance doublet of £ Per (Octo- 
ber 1991). The open circles represent the "first" DAC event that 
reaches a lower vf ymp than its successor (filled diamonds). Al- 
though the two curves reach the same maximum in N co \, they 
are shifted in velocity. 
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Fig. 31. The variation of DAC column density as a function of 
central velocity in the Si IV resonance doublet of 19 Cep. The 
filled diamonds represent the strong DAC event in the August 
1986 dataset; the open circles are the corresponding values for 
the DAC already at vl symp at the beginning of the 1986 cam- 
paign. Clearly, the DAC reaches a maximum in N co i at 0.7 v^. 



4.2. Short-term DAC evolution 

The most striking property of the DAC phenomenon is the 
regularity in recurrence over many years. Time-series analy- 
ses show that the cyclical behaviour of DACs often results in 
a well-determined "wind" period. Although also other periods 
are found (some of them can be identified as harmonics of the 
principal period), the period with strongest power (P w ind) most 
often corresponds to the DAC recurrence timescale t tec . In the 
case of £ Ori t KC (1.6±0.2 day) is also detected in the power 
spectrum, but most power is associated with a period of about 
[6] days. In Table 2 we list -P w i„d derived from the Si IV (orN v) 
resonance lines. The shorter the time span At of the campaign, 
the less reliable the detected periods. When the period is com- 
parable to or longer than At it cannot be trusted. In Table 2 
we marked these periods by putting them in between square 
brackets. 

The DAC recurrence timescale is shorter for stars with 
a higher projected rotational velocity (Henrichs et al 



1988 



Prinja 1988); this conclusion is supported by the results based 



on our larger sample. In Table 3 we compare the recurrence 
timescale of DACs with the estimated maximum rotation pe- 
riod of the star (the inclination of the stellar rotation axis is 
not known). An estimate for the stellar radius is obtained from 
Howarth & Prinja (1989). For all stars t KC is smaller than the 



maximum rotation period of the star. This supports the view 
that stellar rotation determines the characteristic timescale of 
(large-scale) wind variability in O-type stars. Whether t rec is 
equal to the stellar rotation period (or an integer fraction of this) 
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Table 2. Overview of DAC parameters of the studied O stars after modeling the resonance lines of Si IV or N Vindicated with *, 
or both). The dominant period produced by the Fourier analysis is called P w ind- Periods comparable to or longer than the timespan 
At of the dataset are put between square brackets. The listed values are yearly maxima and minima which do not necessarily 
belong to the same DAC event. 



Name 


Yr 


(kms^ 1 ) 


asymp 

(kms -1 ) 


max 
T c 


(kms -1 ) 


(kms" 1 ) 


Armax 
JV col 

(10 13 cm -2 ) 


VNmax 
(km s" 1 ) 


^wind 

(d) 


At 
(d) 






87 
88 
89 
91 


1045±53 
400±141 

1110±71 
660±73 


2438±8 
2439±13 
2141±23 
2319±16 


1.71 ±0.04 
1.09±0.07 
1.05±0.04 
1.23±0.05 


32±6 
54±7 
56±4 
39±11 


1221±110 
2039±202 
691±13 
1265±179 


49.3±1.1 
36.1±0.8 
43.0±0.7 
55.4±0.8 


1617 
1499 
1514 
1771 


1.9±0.3 
1.0±0.1 
1.0±0.1 
2.0±0.2 


3.6 
3.0 
2.5 
4.4 


HD 34656 
slow pattern 


91 


548±25 
1300 


2124±16 
1600 


0.71 ±0.08 
1.07±0.03 


112±16 
52±20 


417±12 
436±12 


26.0±0.8 
24.4±0.6 


1606 
1522 


1.1±0.1 
[4.5±1.6] 


5.1 


AOri A 


92 


1706±21 


2000±15 


0.31 ±0.03 


75±21 


236±20 


4.5±0.3 


1828 


[~4]? 


5.0 


COri A 
N V 


92 


522±19 
551 ±74 


2180±27 
1983±31 


0.41±0.03 
0.21±0.02 


46±12 
31±36 


671 ±62 
454±78 


8.4±0.6 
19.8±2.6 


1129 
1713 


[6.3±2.1] 
[6.1±2.7] 


5.1 


15 Mon* 


91 


1711±38 


2268±65 


0.30±0.08 


43±16 


173±56 


6.7±1.9 


1711 


[5.8±1.8] 


5.3 


68 Cyg 


86 
87 
88 
89 
91 


1141±56 
1200±40 
588±332 
814±41 
840±27 


2499±5 

2448±5 

2565±13 

2419±11 

2488±7 


1.09±0.03 
0.82±0.03 
1.03±0.04 
1.00±0.10 
0.88±0.03 


35±8 
97±15 
26±9 
118±20 
30±10 


562±21 
701 ±24 

1464±414 
998±64 

1258±111 


32.8±0.7 
25.9±0.7 
28.9±1.2 
23.8±0.5 
28.3±0.6 


1587 
1641 
1627 
2187 
1916 


1.4±0.1 

1.3±0.2 

1.3±0.2 

[3.0±0.9] 

1.4±0.2 


6.6 
3.4 
3.1 
2.6 
4.5 


19Cep 
N V 

N V 

N V 


86 
91 
92 


317±94 
545 ±26 
454±63 
1137±181 
771±17 
559±35 


1936±17 
1919±26 
2087±23 
2067±18 
1815±12 
1864±19 


2.28±0.05 
1.13±0.05 
0.78±0.04 
0.80±0.10 
1.14±0.04 
0.90±0.05 


64±18 
42±18 
39±7 
27±19 
99±7 
125±19 


1010±104 
412±100 
846±120 

1185±222 
459±19 
450±69 


68.3±1.1 

118.1±4.7 

29.1±1.0 

67.6±5.2 

34.1±0.8 

101.0±4.7 


1455 
1500 
1494 
1829 
1561 
1647 


[7.0±2.3] 
[7.0±2.3] 
[6.3±2.1] 
[5.6±1.9] 
[4.5±1.9] 
[5.1±1.3] 


6.7 
4.2 
5.6 


ACep 


F91 
091 


703±21 
481 ±22 




1.57±0.07 
1.26±0.10 




484±36 
470±51 






[4.3±1.2] 
1.4±0.2 


5.1 
4.5 


10 Lac* 


92 


710±16 


989±7 


0.38±0.03 


44±10 


166±38 


9.7±0.8 


854 


[6.8±2.3] 


5.0 



will be discussed in section 5. Coordinated Ha and UV obser- 
vations of a few bright O stars have shown that P w ind detected 
in UV resonance lines also appears in the Ha line (Kaper et al. 



1997). On the basis of a large Ha survey Kaper et al. (1998) 



demonstrate that the relation between P w i n( j and the maximum 
rotation period also holds for this larger sample. 

The variability of the inner regions of the stellar wind, 
where Ha is formed, is directly linked to the DAC behaviour in 
UV resonance lines. This is confirmed by the detection of sim- 
ilar variability in subordinate UV lines in £ Per, HD 34656, and 
A Cep. In these stars, the subordinate N IV line at 1718 A with 
P Cygni-type profile (and for A Cep also the He II "Balmera" 
at 1640 A) varies at l ow ve locity in concert with the DACs (see 
also Henrichs et al. 1994 and Paper I). DACs hav e been de- 

As 



1992) 



tected in the N IV line of £ Pup too (Prinja et al. 
these lines have to be formed in regions of higher density, this 
implies that material forming DACs must be present already 
close to the star (within a few i?*). 

Not only the cyclical appearance, but also the change in 
DAC central velocity per unit time is larger for stars with a 
shorter rotation period. This becomes apparent when one com- 
pares the DAC behaviour of the rapid rotators (e.g. £ Per, 
A Cep) with that in the most-likely slowly rotating stars (e.g. 



19 Cep, 10 Lac). This suggests that the "acceleration" of DACs 
is not only a consequence of the line-forming region moving 
away from the star (where the wind speeds are higher); stel- 
lar rotation must play an important role too in explaining the 
acceleration timescale. 

The variation of DAC column density as a function of their 
central velocity (i.e. time) is nicely demonstrated by the appar- 
ently slow rotator 19 Cep and the rapid rotator £ Per in Figs. 
3 1 and 32, respectively. In both cases, the dependence of N co \ 
on v c resembles a parabolic function; a maximum in iV C oi is 
reached half way. As Voo we took the maximum i!c Symp we ob- 
served for the star (Table 1). It turns out that, for a given star, 
not all DACs have the same vl symp . This can vary from cam- 
paign to campaign, or from event to event (e.g. £ Per, Fig. 32). 
Although the absolute velocities are different, the overall shape 
of the variation in N co \ is very similar, even for different stars. 

As said before, the strength of a DAC varies from event 
to event. Sometimes, a weak absorption component disappears 
before reaching the wind terminal velocity. In several stars (e.g. 
£ Per and 68 Cyg) the stronger DACs are accompanied by 
weak absorption components separated by about a quarter of 
the DAC recurrence timescale. In most cases the two compo- 
nents merge when they arrive at the same velocity. 
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Table 3. DAC recurrence timescale compared to the es- 
timated maximum stellar rotation period P max = 50.6 
(vs'mi)^ 1 (R ir / Rq) days for the program stars. An estimate 
for the stellar radius is taken from Howarth & Prinja (1989). 



The stars are ordered according to Pmax- The last column lists 
the proposed stellar rotation period, based on our interpretation 
of the DAC behaviour. 



Name 


i) sin i 




tree 


^max 


Plot 




(kms -1 ) 


(R Q 


(d) 


(d) 


(d) 




68 Cyg 


295 


14 


1.4 


2.4 


2.8 


£Per 


204 


11 


2.0 


2.7 


4.0 


ACep 


214 


19 


1.4 


4.0 


2.8 


HD 34656 


85 


10 


1.1 


6.0 


4.5 


15Mon 


62 


10 


>4.5 


8.2 


>5 


AOri A 


66 


12 


-4? 


9.2 


>4 


19Cep 


90 


18 


~ 5 


10.1 


5 or 10 


COri A 


123 


31 


1.6 


12.8 


6? 


10 Lac 


31 


9 


6.8 


14.7 


7 or 14 





In £ Per and 68 Cyg subsequent (strong) DACs alternate 
in asymptotic velocity. For £ Per the strong components reach 
—2050 or —2300 km s _1 . Since the lifetime of a strong DAC is 
longer than two days (the recurrence timescale), this results in 
"crossings" of successive components; they arrive at the same 
velocity, which does not mean that the material forming the 
two DACs is at the same physical location (though both in the 
line of sight). For 68 Cyg we identified two series of DACs, 
series A reaching an asymptotical velocity of —2450 km s _1 
and series B —2200 km s _1 . A difference with £ Per is that the 
time interval between series B and A is 0.5 d, which is much 
shorter than the recurrence timescale of 1 .4 d (the distance in 
time between subsequent A's or B's, as well as pairs B A). 
The difference between these sequences of DACs is visualized 
as follows: 





series A 




^asymp 


series B — * 


v 2 

asymp 


£Per 


A B 
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B A 


BAB 


A ... 


68 Cyg 
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B A 


B ... 


HD 34656 


A B 
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B A 


B B B 
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time — > 



The above mentioned DAC properties make up for the ob- 
served "pattern" of variability that is characteristic for the star. 
Sometimes it seems to be more complicated; e.g. in HD 34656 
two different patterns can be recognized, a fast one and a slow 
one. More extended time series, like those obtained during the 
IUE MEGA campaign (Massa et al. |1995|), also indicate that 



the regular variability pattern is modulated by variations oc- 
curing over a longer, but eventually related timescale. A good 
example is provided by the MEGA campaign data on £ Pup 



(Howarth et al. 1995). These authors interpret the longer (5.2- 
day) period as the stellar rotation period; the shorter (19-hour) 
recurrence timescale of the DACs would, however, not be an in- 
teger fraction (6.5±0.9) of the rotation period. If the latter were 



true for this star, stellar rotation alone would not be sufficient 
to explain the cyclical recurrence of DACs. 

4.3. Edge variability 

In saturated P Cygni profiles wind variability manifests itself in 
significant shifts in velocity of the steep blue edge (cf. Paper I). 
In section 3 it is shown that close to the high-velocity edges of 
the UV resonance lines often a longer period is detected than 
the period exhibited by the DACs. We propose that the edge 
variability is directly related to the DAC behaviour. The obser- 
vations suggest that the slow modulation of the high-velocity 
edges is the result of the fact that only the DACs with the high- 
est asymptotic velocity have an impact on the position of the 
edge. The case of £ Per provides the strongest support to this 
interpretation: the edge variability has a period of about 4 days 
(i.e., twice the DAC recurrence timescale), while only one out 
of two DACs reaches an asymptotic velocity exceeding —2050 
km s — 1 . 

The combination of a slow and a more rapid variability pat- 
tern (such as observed in HD 34656and ( Pup) and the subse- 
quent modulation in strength of the DACs might give rise to 
edge variability in saturated profiles on the timescale of the 
slow pattern. Therefore, only extended time series (such as 
those obtained in the IUE MEGA campaign, cf. Massa et al. 



1995) might enable confirmation of this conjecture. 



4.4. Phase diagrams 

The phase diagrams included in this paper show the variation 
of the sinusoid's phase at the dominant wind variability period 
with position in the line. The observed decrease in phase to- 
wards higher velocities is due to the acceleration of the DACs 
through the line profile. This supports our view that P w i n( j is 
equal to the DAC recurrence timescale. The large width of the 
DACs when they appear at low or intermediate velocity would 
explain the relatively small change in phase observed at these 
velocities. 

Owocki et al. (1995) discuss the phase "bowing" discov- 
ered in the phase diagrams belonging to P w ind detected in the 



B0.5 lb star HD 64760 (Prinja et al. |1995| , Fullerton et al. |1997[ ). 
Although, according to these authors, in this star P w i n d (1.2 or 
2.4 d) is caused by a modulation in flux at low and intermediate 
velocities rather than by DACs, the decrease of phase towards 
both higher and lower velocities (which appears as a bow in the 
phase diagram) is a strong diagnostic of curved wind structures 
(like co-rotating interacting regions, section 5). DACs are de- 
tected at higher velocities in the resonance lines of HD 64760 
and are variable on a timescale longer than P w i nc |. 

The phase diagrams for the 2-day period detected in the 
£ Per October 1991 dataset clearly show evidence for phase 
bowing. A maximum in phase occurs around —1200 km s -1 , 
i.e. about the central velocity at which the DACs usually ap- 
pear (Fig. 8). The maximum phase lag A(j> is about one radian, 
similar to the phase lag observed in HD 64760. The increasing 
strength of the DAC might explain the subsequent decrease in 
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phase towards higher and lower velocities. The further acceler- 
ation and narrowing of the DACs is reflected by the increasing 
slope of the phase towards higher velocities. This description 
of the phase behaviour is only qualitative; detailed modeling of 
the phase diagram falls beyond the scope of this paper. An im- 
portant difference with HD 64760 is that the DACs themselves 
seem to be responsible for the observed phase bowing. 

Phase bowing is not evident in the phase diagrams for the 
other O stars in our sample. It has only been found in £ Per, 
which is perhaps not surprising, because the other stars in our 
sample do not exhibit such strong wind variability down to the 
lowest wind velocities. 

4.5. Flux modulation and DAC behaviour 

The "classical" picture we have of DAC behaviour, and which 
is supported by the observations analysed in this paper, is that 
broad absorption components appear at low velocity which 
evolve towards higher velocity while narrowing. The applica- 
tion of Fourier analyses resulted in the determination of well- 
constrained periods, not only at low velocity, but over the full 
width of the lines. To check the outcome of the Fourier analy- 
sis we constructed diagrams showing the variation of flux as 
a function of time in a each velocity bin. It indeed appears 
that the flux varies sinusoidally with time, especially at low 
and intermediate velocities. In some cases the Fourier analysis 
invokes an additional sinusoid to accommodate the weak pre- 
cursor of the next DAC (resulting in the addition of a shorter- 
period sinusoid) or a gradual change (longer period). The latter 
occurs for example when two subsequent DACs reach a differ- 
ent asymptotic velocity. 

The analysis of two spectroscopic time series of UV reso- 
nance lines of the BO. 5 lb supergiant HD 64760 lead Fullerton 
et al. (1997) to conclude that in this star a distinction can be 



made between a sinusoidal modulation in flux at low and in- 
termediate wind velocities and variations due to the propaga- 
tion of DACs at higher velocities. A 2.4-day (or 1.2-day) peri- 
odicity is connected with the modulation which clearly shows 
the occurrence of phase bowing; the DACs vary on a longer 
timescale. Thus, in HD 64760 two different kinds of wind vari- 
ability are observed: flux modulation and DAC propagation. 



Because of the observed phase bowing, Fullerton et al. (1997) 



interpret the modulation to be due to CIRs, whereas the origin 
of the DACs is not clear. 

Our results on £ Per indicate a different picture: in this star 
the flux modulation and phase bowing are caused by the DACs 
themselves. That the propagation of DACs leads to a sinusoidal 
modulation of the flux in a given velocity bin is due to the shape 
of the DACs, which can be accurately fit with an exponential 
Gaussian. At low velocity these gaussians have widths of sev- 
eral hundreds km/s. At high velocity they become very narrow. 
When the DAC moves through the resonance line, the flux in 
a given velocity bin changes gradually, following the Gaussian 
shape of the DAC. At high velocities, where the DACs are nar- 
row, the flux changes more irregularly, because the DACs "pile 
up" at about the same asymptotic velocity before they vanish. 



This might explain why the flux gradually changes when a 
DAC moves through a velocity bin, but not why this modula- 
tion is sinusoidal. In our interpretation, this is because not only 
the recurrence timescale, but also the acceleration timescale of 
the DAC is proportional to the rotation period of the star. Since 
these two timescales are identical, a sinusoidal shape of the flux 
versus time plots is expected. Thus, the DAC behaviour in O- 
star winds can result in a sinusoidal modulation of the flux in 
each velocity bin. 

With our definition of the template, the maximum of the 
sinusoidal modulation corresponds to the "undisturbed, under- 
lying wind profile", so that all variations are by definition due 
to a dditio nal absorption. From the point of view of Fullerton et 
al. ( |1997| ), in the case of HD 64760 the template might be bet- 
ter defined as the zero-point of the modulation, i.e. the "average 
wind profile". With such a template, it is not possible to model 
the DACs they observe in the resonance lines of HD 64760 at 
high velocity as absorption components, which is the reason 
for our approach. 

4.6. Long-term variability 

The O stars £ Per, 68 Cyg, 19 Cep, and A Cep were observed 
several times in the period 1986-1992. The variability pattern, 
characteristic for a given star, can be recognized from year to 
year, though detailed changes occur. The same wind variabil- 
ity timescale P w i nc j is detected every year. Sometimes a differ- 
ent period gives the highest peak in the power spectrum (£ Per 
1988 and 1989, 68 Cyg 1989), but this might be due to the 
short time coverage of these datasets. For A Cep there seems to 
be a clear difference between the observed variability timescale 
in October 1989, February 1991, and October 1991. In 1989 a 
period of [2.5±0.8 d] is dominant, in February 1991 a period 
of [4.3±1.2 d] (as well as 1.4±0.1 d), and in October 1991 a 
period of 1.4±0.2 d. If stellar rotation is setting the timescale 
of wind variability, one would not expect this timescale to vary 
from year to year. What could in principle vary is the number 
of features observed per rotation period. In the case of A Cep 
the difference in observed periods could be explained as dif- 
ferent integer fractions of the stellar rotation period: 1.4, 2.8, 
4.2 d; or 1.3, 2.6, 3.9 d. In section 5 we will address the point 
whether the DAC recurrence timescale can be used to deter- 
mine the stellar rotation period. 

Although the characteristic variability timescale remains 
the same, the DAC strengths and asymptotic velocities exhibit 
changes from campaign to campaign. The DAC behaviour in 
68 Cyg seems to be the most constant over the years; clear dif- 
ferences are observed for £ Per (Fig. 6), 19 Cep, and A Cep. 
This indicates that the underlying clock might be the same, 
but that at least for some stars in our sample the pattern of 
large-scale structure in the stellar wind varies on a timescale 
of months to years. 
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5. Towards an interpretation 

Several models have been put forward to explain the observed 
properties of DACs. The expan sion o f a high-density layer in 
the stellar wind (Henrichs et al. 1983 ) gives a simple explana- 
tion for the time-dependent behaviour of DACs, but it gives no 
clue to the origin of the high-density layer. Underhill & Fahey 
( 1984 ) tried to explain the existence and behaviour of DACs 
by postulating the existence of "closed" and "open" magnetic 
loops above the stellar surface from which parcels of gas are 
released in addition to a uniformly emitted steady stellar wind. 
A nice property of this model is that DACs are always expected 
to appear at the same rotational phase as long as the magnetic 
footpoints remain at the same location. 

DACs are found at the same velocity in reso nance lines 
formed by different ions (e.g. Lamers et al. 1982); therefore, 



DACs most-likely originate from regions which have a higher 
density (or smaller velocity gradient) than the ambient wind. 
These regions have to be geometrically extended, because: (1) 
practically all O stars show DACs; (2) they have to cover a sig- 
nificant fraction of the stellar disk to be observable. 

Stellar rotation plays a dominant role in setting the 
timescale of variability; a difference in flow time scale 
(R+/voo) cannot explain the observed dependence on vsini. 
Not only the cyclical recurrence, but also the acceleration of 
DACs is faster for more rapidly rotating stars. This suggests 
that the DAC forming region rotates into and out of the line of 
sight. The observed increase and subsequent decrease of A^oi 
supports this picture (Figs. 31 and 32). 

5.1. Corotating Interacting Regions and DACs 

Most-promising is the Corotating Interacting Regions model 
proposed by Mullan (1984, 1986), in analogy with the coro- 
tating interacting regions in the solar wind which result from 
the interaction of fast and slow wind streams. Recently, this 
model has been worked out by Cranmer & Owocki (1990 
who performed two-dimensional hydrodynamical simulations 
of corotating stream structure in the wind of a rotating O star. 
The corresponding spiral-shaped, large-scale structure through 
which the wind material is flowing, was proposed by Prinja 



& Howarth (1988) based on observations of 68 Cyg. Simulta- 



neous ultraviolet and Ha spectroscopy of a number of bright 
O stars included in this study strongly support the CIR model 
(Kaper et al. 1997 see also Fullerton et al. 1997). Cranmer 



& Owocki (1996) conclude from their model that the DACs 



are not formed in the CIRs themselves, but originate in the so- 
called radiative-acoustic kinks trailing the CIRs. The strong ve- 
locity gradient related to this kink results in a relativel y large 
contribution to the Sobolev optical depth. Kaper et al. ( 1997 ) 
show that a small phase lag might be present between the wind 
variations in the Ha line with respect to those observed in the 
UV resonance lines; this would be consistent with the Ha vari- 
ations originating in the CIR and the DACs in the trailing kink. 

Phase bowing suggests that the structures are curved, ex- 
plaining why material at a given radial velocity rotates out of 



or into the line of sight first, followed by material moving at 
higher and lower velocities (Owocki et al. 1995, Cranmer & 
Owocki 1996| ). All nine O stars with DACs show a decrease in 
phase towards higher velocities; only one (£ Per) also shows 
a lag in phase towards lower velocities. The location in the 
wind where the fast flow first meets the slow stream depends on 
the equatorial rotation speed v rot and the difference in (termi- 



nal) velocity between the fast and slow stream (Mullan 1984) 
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Inwards of rj, the stellar wind might be rather smooth, so that 
variability is not detected at low wind velocities. Also, the 
"contrast" between the fast and slow stream (angular separa- 



tion Atp and v 



fast 

oo 



v ) is likely to vary from star to star. For 



large it might well be that only a decrease in phase towards 
higher velocities will be detected. Furthermore, the aspect an- 
gle might be an important factor; e.g. the decrease in phase to- 
wards lower velocities is only observed when the line of sight 
includes the equatorial plane. 

The CIR model only works for rotating stars with emerging 
wind flows of different kinematic structure. Comparison be- 
tween model and observations suggest that wind streams with 
relatively slow speed interact with the faster ambient wind, 
in stead of hi gh-speed streams in a slow wind (Cranmer & 
Owocki 1996 ). Therefore, the asymptotic velocity reached by 
DACs might be somewhat lower than the terminal velocity of 
the ambient wind. In some stars we measure systematically dif- 
ferent values of vf ymp from event to event. This might be due 
to a difference in aspect angle, contrast, or both. 

5.2. The origin of CIRs: non-radial pulsations or surface 
magnetic fields ? 

In order to work, this model needs a certain structure imposed 
at the stellar surface to produce flows with different kinemati- 
cal properties. Two physical mechanisms that in principle could 
cause fast and slow streams in the wind are: (i) non-radi al pul - 
sations (NRP), or (ii) surface magnetic fields. Henrichs (1984) 
pointed out that the occurrence of DACs might be related to 
the presence of NRP at the stellar surface. NRP are detected 
(or suspected) in several O stars (cf. Baade 1988, Fullerton et 
al. 1996). The associated temperature and velocity fields (on 



the order of the sound speed) could be responsible for the mod- 
ulation of the boundary conditions at the base of the radiation- 
driven wind. 

Our working hypothesis is that the outflow properties at the 
base of the wind are influenced by a magnetic field anchored 
in the star, which results in regions with different outflow prop- 
erties rotating with the star. The presence of magnetic fields 
in early-type stars is difficult to prove with direct observa- 
tional methods. Magnetic fields with a strength of less than 100 
Gauss are theoretically predicted (Maheswaran & Cassinelli 



1992), but not (y et) observable (Landstreet 1992). Recently, 



Stahl et al. (|1996j) argued that a surface magnetic field might be 
present in the 07 V star 9 1 Ori C; phase-locked photospheric 
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and stellar-wind variations strongly suggest that this star is an 
oblique magnetic rotator. An upper limit of 70 G for the (lon- 
gitudinal) magnetic field strength of £ Per has been derived by 



Henrichs et al. ( 1998 ). Recent attempts to measure the longitu- 



dinal magnetic field strength in a sample of 10 fairly slowly ro- 



tating O stars by Mathys (1998) and collaborators did not result 



in a definite detection (la error bars of the order of 250-300 G). 
Some marginal measurements (between 2 and 3a) suggest that 
weak fields might nevertheless be present, which could become 
detectable through a fairly modest improvement of the mea- 
surement accuracies. 

Early-type stars are relatively short-lived objects and born 
out of interstellar matter that contains a magnetic field. The star 
with the strongest magnetic field known, "Babcock's star" HD 
215441, is a B3 V star with a surface dipole field strength of 
34 kGauss (Babcock I960, Preston 1971). This demonstrates 



that early-type stars can have very strong magnetic fields (see 
also Bohlender 1994). These fields are expected to be frozen-in 
into the radiative envelopes of the stars (see, for example Mes- 



tel 1975 1. Furthermore, the end-products of early-type stars 



are mostly neutron stars with a strong magnetic field, which 
strength can be predicted by assuming magnetic flux conserva- 
tion during the collapse of the core and the presence of a sur- 
face magnetic field of strength 10-100 G in the progenitor. Al- 
ternative scenarios for the formation of neutron-star magnetic 
fields exist as well. 



5.3. Can t rec be used to measure P m t? 

When we adopt the CIR model in combination with an an- 
chored surface magnetic field to explain the presence and cycli- 
cal behaviour of DACs, the characteristic timescale of variabil- 
ity in the winds of O-type stars is the rotation period of the star. 
Our observations suggest that in several cases (e.g. £ Per and 
68 Cyg where successive DACs have different asymptotical ve- 
locities) the rotation period of the star is more likely equal to 
twice (or an integral number times) the recurrence timescale 
of DACs. This means that DACs are not only a diagnostic for 
deriving the terminal velocity of a stellar wind, but they can 
in principle be used to derive the rotation period of the star. In 
Table 3 we list, based on the interpretation of the cyclical DAC 
behaviour, the best estimate for the rotation period of the star. 
Given an estimate of the stellar radius, the inclination of the 
rotation axis can thus be derived from v sin i. 

In two cases (£ Per and 68 Cyg) the proposed rotation 
period exceeds the estimated maximum rotation period. It is 
likely that these rapidly rotating stars (on the basis of the mea- 
sured v sin i) have a sin i value close to unity, and the radius 
of these stars could be slightly underestimated. For £ Per a ra- 
dius of 15i?0 instead of lli? Q would increase the maximum 



rotation period to 4 days (see Puis et al. 1996). Walborn et al 



(1985) note that £ Per is of luminosity class II instead of III 



which supports this conclusion. The uncertainty in the spectral 
class of 68 Cyg (denoted by the colon) also indicates that the 
estimated radius might not be accurate. 



5.4. Summary and future work 

The DAC behaviour in the UV resonance lines of a sample of 
10 bright O stars has been analysed in detail. The cyclical ap- 
pearance and subsequent evolution of DACs is interpreted in 
terms of a model invoking fast and slow streams which inter- 
act due to the rotation of the underlying star. The interacting 
regions (CIRs) corotate with the star while the wind material is 
flowing through them. The observations suggest that the stellar 
wind includes more than one CIR, most likely two (cf. Kaper et 
al. 1997). The curvature of the CIRs can in principle explain the 
observed phase behaviour (bowing), which is detected in one of 
the O stars in our sample. The wind variability periods we de- 
rive are a direct measure of the DAC recurrence timescales and, 
in our interpretation, equal to an integer fraction (most likely 
1/2) of the stellar rotation period. 

The key issue is the origin of the CIRs: are non-radial pul- 
sations, surface magnetic fields, or other physical mechanisms 
responsible for the surface structure creating fast and slow wind 
streams? This is the subject of our current investigations. 
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Fig. 6. £ Per 07.5 III(n)((f)): An overview of the DAC behaviour in the Si IV resonance lines over a period of 4 years. The 
quotient spectra (middle panel) are shown as a "dynamic" spectrum in the grey-scale panel below. Outstanding is the similarity 
of the DAC patterns, though detailed changes do occur. The significance of variability (Paper I) is indicated in the upper panels 
(thick line), together with the mean spectrum of that particular year (thin line). 
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Fig. 7. £ Per September 1987 and October 1988: Parameters of DACs in the Si IV resonance doublet. The central velocity v c , 
central optical depth r c , and width v t of the components are plotted as a function of time (in days). The figures belonging to 
datasets of the same star are plotted on identical scale. The column density iV co i is calculated with Eq. (10) and expressed in 
units of 10 13 cm~ 2 . In the top panel the total equivalent width (EW) of the Si IV quotient spectra is given in km s _1 . The points 
suggested to belong to the same DAC event are interconnected. 
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Fig. 8. £ Per October 1989 and October 1991 (as Fig. 7): DACs in the Si IV resonance doublet. 
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HD 34656 Si IV DAC fits February 1991 
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Fig. 13. As Fig. 7: parameters of DACs in the Si IV doublet at 1400 A of HD 34656 07 11(f) in February 1991. The DAC 
behaviour can be interpreted in terms of a "slow" and a "fast" pattern, the latter with a recurrence timescale close to one day. 
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Fig. 16. As Fig. 7: ( Ori 09.7 lb November 1992: fit parameters of DACs in the N V (left panel) and Si IV profile (right panel). 
A new DAC appears about every 1.5 days. 




Fig. 17. As Fig. 7: 15 Mon 07 V((f)) February 1991: two DACs are detected in the N V time series. The DAC appearing at day 9 
has a width v t < 150 km s _1 ; this narrow width has also been observed for the DACs in A Ori and 10 Lac. 
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68 Cyg Si IV DAC fits 
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Fig. 18. 68 Cyg 07.5 I II:n((f )) August 1986 (as Fig. 7): parameters of DACs in the Si IV doublet. This dataset was also analyzed 
by Prinja & Howarth (1988). During the first three days time coverage is very sparse, so that the different DACs could not be 
labeled. Close inspection of the available datasets (cf. Figs. 19 and 20) results in the recognition of typical DAC patterns. On this 
basis the DACs are labeled with symbols to highlight the pattern. 
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68 Cyg Si IV DAC fits October 1989 



October 1991 
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Fig. 20. 68 Cyg October 1989 and October 1991 (as Fig. 7): DAC parameters in the Si IV doublet. The October 1991 dataset has 
the longest time coverage (4.5 d) and includes 40 spectra. 
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19 Cep August 1986 Si IV DAC fits N V DAC fits 
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Fig. 23. 19 Cep 09.5 lb August 1986 (as Fig. 7): DAC model parameters for the Si IV (left) and nv (right) doublets. The two 
profiles give consistent results and show slowly evolving DACs. The DAC appearing at day 7 is the strongest component we 
encountered in our O star dataset. 
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Fig. 26. A Cep 06 I(n)fp in October 1991 (as Fig. 7): DAC parameters resulting from profile fitting of the Si IV doublet. The 
velocity domain from — 1 600 to —2000 km s ~ 1 is disturbed by large fluctuations, because of saturation of the absorption profile. 
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Fig. 27. A Cep 06 I(n)fp in October 1991: Grey-scale representations showing the time evolution of the Si IV resonance doublet 
and the He II and N IV subordinate lines. A template spectrum (upper panels) was used to produce residual spectra showing the 
development and acceleration of 4 DAC events with better contrast. The variability amplitude is illustrated by a thick line in the 
upper panel. The first signs of DACs are observed at a velocity of —200 km s _1 in these subordinate lines. 
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